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1-^^ ' Abstract 



We consider the physics potential of CERN based neutrino oscihation experiments 
consisting of a Beta Beam (/3B) and a Super Beam (SPL) sending neutrinos to MEM- 
5^ \ PHYS, a 440 kt water Cerenkov detector at Frejus, at a distance of 130 km from 

CERN. The ^13 discovery reach and the sensitivity to CP violation are investigated, 
including a detailed discussion of parameter degeneracies and systematical errors. For 
SPL sensitivities similar to the ones of the phase II of the T2K experiment (T2HK) 
are obtained, whereas the (3^ may reach significantly better sensitivities, depending 
on the achieved number of total ion decays. The results for the CERN-MEMPHYS 
experiments are less affected by systematical uncertainties than T2HK. We point out 
that by a combination of data from /3B and SPL a measurement with antineutrinos is 
not necessary and hence the same physics results can be obtained within about half of 
the measurement time compared to one single experiment. Furthermore, it is shown 
how including data from atmospheric neutrinos in the MEMPHYS detector allows to 
resolve parameter degeneracies and, in particular, provides sensitivity to the neutrino 
mass hierarchy and the octant of ^23- 
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1 Introduction 



In recent years strong evidence for neutrino oscillations has been obtained in solar [1], at- 
mospheric [2,3], reactor [4], and accelerator [5] neutrino experiments. The very near fu- 
ture of long-baseline (LBL) neutrino experiments is devoted to the study of the oscillation 
mechanism in the range of Am^^ ~ 2.4 x 10~^ eV^ indicated by atmospheric neutrinos using 
conventional beams. Similar as in the K2K experiment in Japan [5], the presently running 
MINOS experiment in the USA [6] uses a low energy beam to measure Amg^ by observing the 
~^ disappearance probability, while the OPERA [7] experiment will be able to detect 
Pr appearance within the high energy CERN-Gran Sasso beam [8]. If we do not consider the 
LSND anomaly [9] that will be further studied soon by the MiniBooNE experiment [10], all 
data can be accommodated within the three flavor scenario (see Refs. [11,12] for recent global 
analyses), and neutrino oscillations are described by two neutrino mass-squared differences 
(Am^i and Arrig]^) and the 3x3 unitary Pontecorvo-Maki-Nakagawa-Sakata (PMNS) lepton 
mixing matrix [13] with three angles (^12,^13,^23) and one Dirac CP phase Sep- 

Future tasks of neutrino physics are an improved sensitivity to the last unknown mixing 
angle, ^^13, to explore the CP violation mechanism in the leptonic sector, and to determine 
the sign of Amli which describes the type of the neutrino mass hierarchy (normal, Am§]^ > 
or inverted, Amg^^ < 0). The present upper bound on ^13 is dominated by the constraint 
from the Chooz reactor experiment [14]. A global analysis of all data yields sin^ 2^13 < 0.082 
at 90% CL [12]. A main purpose of upcoming reactor and accelerator experiments is to 
improve this bound or to reveal a finite value of ^13. In reactor experiments, one uses z/g 
in disappearance mode and the sensitivity is increased with respect to present experiments 
by the use of a near detector close to the reactor [15]. In accelerator experiments, the first 
generation of so-called Super Beams with sub-mega watt proton drivers such as T2K (phase- 
I) [16] and NOz/A [17], the appearance channel v^. is explored. This next generation of 

reactor and Super Beam experiments will reach sensitivities of the order of sin^ 2^13 < 0.01 
(90% CL) within a time scale of several years [18]. Beyond this medium term program, there 
are several projects on how to enter the high precision age in neutrino oscillations and to 
attack the ultimate goals hke the discovery of leptonic CP violation or the determination of 
the neutrino mass hierarchy. In accelerator experiments, one can extend the Super Beam 
concept by moving to multi-mega watt proton drivers [16, 19-21] or apply novel technologies, 
such as neutrino beams from decaying ions (so-called Beta Beams) [22,23] or from decaying 
muons (so-called Neutrino Factories) [23,24]. 

In this work we focus on possible future neutrino oscillation facilities hosted at CERN, 
namely a multi-mega watt Super Beam experiment based on a Super Proton Linac (SPL) [25] 

and a7 = 100 Beta Beam (/3B) [26]. These experiments will search for and ''v e-^'^v jx 

appearance, respectively, by sending the neutrinos to a mega ton scale water Cerenkov 
detector (MEMPHYS) [27], located at a distance of 130 km from CERN under the Frejus 
mountain. Similar detectors are under consideration also in the US (UNO [28]) and in Japan 
(Hyper-K [16,29]). We perform a detailed analysis of the SPL and /3B physics potential, 
discussing the discovery reach for ^13 and leptonic CP violation. In addition we consider 
the possibility to resolve parameter degeneracies in the LBL data by using the atmospheric 
neutrinos available in the mega ton detector [30]. This leads to a sensitivity to the neutrino 
mass hierarchy of the CERN-MEMPHYS experiments, despite the rather short baseline. 
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SPL 



T2HK 



Detector mass 
Baseline 



Running time (y + v) 
Beam intensity 



440 kt 440 kt 

130 km 130 km 

5 + 5 yr 2 + 8 yr 



440 kt 
295 km 
2 + 8 yr 
4 MW 

2% 

2% 



5.8 (2.2) • lO^s He (Ne) dcys/yr 4 MW 



Systematics on signal 
Systematics on backgr 



2% 2% 
2% 2% 



Table 1: Summary of default parameters used for the simulation of the /3B, SPL, and T2HK experiments. 

The physics performances of /SB and SPL are compared to the ones obtainable at the second 
phase of the T2K experiment in Japan, which is based on an upgraded version of the original 
T2K beam and the Hyper-K detector (T2HK) [16]. 

The outline of the paper is as follows. In Sec. [2] we summarize the main characteristics 
of the /3B, SPL, and T2HK experiments and give general details of the physics analysis 
methods, whereas in Sec. [3] we describe in some detail the MEMPHYS detector, the /5B, the 
SPL Super Beam, and our atmospheric neutrino analysis. In Sec. |l]we review the problem of 
parameter degeneracies and discuss its implications for the experiments under consideration. 
In Sec. Owe present the sensitivities to the "atmospheric parameters" ^23 and Am|^, the 
^13 discovery potential, and the sensitivity to CP violation. We also investigate in some 
detail the impact of systematical errors. In Sec. we discuss synergies which are offered 
by the CERN-MEMPHYS facilities. We point out advantages of the case when /3B and 
SPL are available simultaneously, and we consider the use of atmospheric neutrino data 
in MEMPHYS in combination with the LBL experiments. Our results are summarized in 



2 Experiments overview and analysis methods 

In this section we give the most important experimental parameters which we adopt for the 
simulation of the CERN-MEMPHYS experiments /?B and SPL, as well as for the T2HK 
experiment in Japan. These parameters are summarized in Tab. [H For all experiments the 
detector mass is 440 kt, and the running time is 10 years, with a division in neutrino and 
antineutrino running time in such a way that roughly an equal number of events is obtained. 
We always use the total available information from appearance as well as disappearance 
channels including the energy spectrum. For all three experiments we adopt rather optimistic 
values for the systematical uncertainties of 2% as default values, but we also consider the case 
when systematics are increased to 5%. These errors are uncorrelated between the various 
signal channels (neutrinos and antineutrinos) , and between signals and backgrounds. 

A more detailed description of the CERN-MEMPHYS experiments is given in Sec. El 
For the T2HK simulation we use the setup provided by GLoBES [31] based on Ref. [32], 
which follows closely the LOI [16]. In order to allow a fair comparison we introduce the 
following changes with respect to the configuration used in Ref. [32]: The fiducial mass is 
set to 440 kt, the systematical errors on the background and on the z/g and z/g appearance 
signals is set to 2%, and we use a total running time of 10 years, divided into 2 years of data 
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Table 2: Number of events for appearance and disappearance signals and backgrounds for the /3B, SPL, 
and T2HK experiments as defined in Tab. [T] For the appearance signals the event numbers are given for 
several values of sin^ 2^13 and (5cp and 7r/2. The background as well as the disappearance event numbers 
correspond to 6*13 = 0. For the other oscillation parameters the values of Eq. Ilj are used. 



taking with neutrinos and 8 years with antineutrinos. We include an additional background 
from the v^J^ — >■ z/g (f^ u^) channel in the neutrino (antineutrino) mode. Furthermore, we 
use the same CC detection cross section as for the /5B/SPL analysis [33]. For more details 
see Refs. [16,32]. 

Recently the idea was put forward to observe the T2K beam with a second detector placed 
in Korea [34-36]. In the so-called T2KK setup the Hyper-Kamiokande detector is split into 
two detectors of 270 kt each, one of them is located in Korea at a distance of about 1050 km 
from the source, and the other is placed at Kamioka at a distance of 295 km. Here we will 
confine ourselves to the standard T2HK setup, since the purpose of our work is not a T2HK 
optimization study investigating various configurations for that experiment. In contrast, here 
T2HK mainly serves as a point of reference to which we compare the CERN-MEMPHYS 
experiments. For this aim we prefer to stick to the "minimal" one-detector configuration at a 
relatively short baseline, since two-detector setups with very long baselines clearly represent a 
different class of experiments whose consideration goes beyond the scope of the present work. 
Nevertheless we will comment briefly also on T2KK performances obtained in Refs. [34,35]. 

In Tab. [2] we give the number of signal and background events for the experiment setups 

as defined in Tab. [H For the appearance channels for the /5B and ^z^'^— >'z^'e for 

SPL and T2HK) we give the signal events for various values of 6*13 and 6cp- The "signal" 
events for ^13 = are appearance events induced by the oscillations with Amji. The value 
sin^ 26*13 = 10"'^ corresponds roughly to the sensitivity limit for the considered experiments, 
whereas sin^ 26*13 = 10"^ gives a good sensitivity to CP violation. This can be appreciated 
by comparing the values of z/ and u appearance events for 6cp = and Tr/2. In the table the 
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background to the appearance signal is given for ^13 = 0. Note that in general the number 
of background events depends also on the oscillation parameters, since also the background 
neutrinos in the beam oscillate. This effect is consistently taken into account in the analysis, 
however, for the parameter values in the table the change in the background events due to 
oscillations is only of the order of a few events. 

The physics analysis is performed with the GLoBES open source software [31], which 
provides a convenient tool to simulate long-baseline experiments and compare different fa- 
cilities in a unified framework. The experiment definition (AEDL) files for the (3B and SPL 
simulation with GLoBES are available at Ref. [31]. In the analysis parameter degeneracies 
and correlations are fully taken into account and in general all oscillation parameters are 
varied in the fit. To simulate the "data" we adopt the following set of "true values" for the 
oscillation parameters: 

Aml^ = +2.4 X 10-3 eV^ , sin^ 623 = 0.5 , . . 

Aml^ = 7.9 X 10-5 eV^ , sin^ ^12 = 0.3 , ^ ^ 

and we include a prior knowledge of these values with an accuracy of 10% for 612, 623, Arrig^, 
and 4% for Amli at Icr. These values and accuracies are motivated by recent global fits to 
neutrino oscillation data [11,12], and they are always used except where explicitly stated 
otherwise. 



3 The CERN-MEMPHYS experiments 
3.1 The MEMPHYS detector 

MEMPHYS (MEgaton Mass PHYSics) [27] is a mega ton class water Cerenkov detector in the 
straight extrapolation of Super-Kamiokande, located at Frejus, at a distance of 130 km from 
CERN. It is an alternative design of the UNO [28] and Hyper-Kamiokande [29] detectors and 
shares the same physics case, both from the non-accelerator domain (nucleon decay, super 
nova neutrino detection, solar neutrinos, atmospheric neutrinos) and from the accelerator 
domain which is the subject of this paper. A recent civil engineering pre-study to envisage 
the possibly of large cavity excavation located under the Frejus mountain (4800 m.e.w.) near 
the present Modane underground laboratory has been undertaken. The main result of this 
pre-study is that MEMPHYS may be built with present techniques as a modular detector 
consisting of several shafts, each with 65 m in diameter, 65 m in height for the total water 
containment. A schematic view of the layout is shown in Fig. [TJ For the present study we 
have chosen a fiducial mass of 440 kt which means 3 shafts and an inner detector of 57 m 
in diameter and 57 m in height. Each inner detector may be equipped with photo detectors 
(81000 per shaft) with a 30% geometrical coverage and the same photo-statistics as Super- 
Kamiokande (with a 40% coverage). In principle up to 5 shafts are possible, corresponding 
to a fiducial mass of 730 kt. The Frejus site offers a natural protection against cosmic rays 
by a factor 10^. If not mentioned otherwise, the event selection and particle identification 
are the Super-Kamiokande algorithms results. 
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3.2 The 7 = 100 X 100 baseline Beta Beam 



The concept of a Beta Beam (/3B) has been introduced by P. Zucchelli in Ref. [22]. Neutrinos 
are produced by the decay of radioactive isotopes which are stored in a decay ring. An 
important parameter is the relativistic gamma factor of the ions, which determines the energy 
of the emitted neutrinos. jSB performances have been computed previously for 7(^He) = 
66 [26], 100 [37-39], 150 [39], 200 [40], 350 [39], 500 [40,41], 1000 [40], 2000 [41], 2488 [42]. 
Reviews can be found in Ref. [43], the physics potential of a very low gamma j3B has been 
studied in Ref. [44]. Performances of a /?B with 7 > 150 are extremely promising, however, 
they are neither based on an existing accelerator complex nor on detailed calculations of 
the ion decay rates. For a CERN based (3B, fluxes have been estimated in Ref. [45] and a 
design study is in progress for the facility [46]. In this work we assume an integrated flux 
of neutrinos in 10 years corresponding to 2.9 ■ 10^^ useful ^He decays and 1.1 ■ 10^^ useful 
^^Ne decays. These fluxes have been assumed in all the physics papers quoted above, and 
they are two times higher than the baseline fluxes computed in Ref. [45]. These latter fluxes 
suffer for the known limitations of the PS and SPS synchrotrons at CERN, ways to improve 
them have been delineated in Ref. [47]. 

The infrastructure available at CERN as well as the MEMPHYS location at a distance 
of 130 km suggest a 7-factor of about 100. Such a value implies a mean neutrino energy of 
400 MeV, which leads to the oscillation maximum at about 200 km for Am|]^ = 2.4x 10~^ eV^. 
We have checked that the performance at the somewhat shorter baseline of 130 km is rather 
similar to the one at the oscillation maximum. Moreover, the purpose of this paper is to 
estimate the physics potential for a realistic set-up and not to study the optimization of 
the j3B regardless of any logistic consideration (see, e.g., Refs. [39,40] for such optimization 
studies) . 

The signal events from the —>■ neutrino and antineutrino appearance channels in 
the (3B are charged current (CC) events. The Nuance v3r503 Monte Carlo code [33] is 
used to generate signal events. The selection for these events is based on standard Super- 
Kamiokande particle identification algorithms. The muon identification is reinforced by 
asking for the detection of the Michel decay electron. The neutrino energy is reconstructed 
by smearing momentum and direction of the charged lepton with the Super-Kamiokande res- 
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Figure 2: Energy resolution for interactions in the 200-300 MeV energy range. The quantity displayed 
is the difference between the reconstructed and the true neutrino energy. 



olution functions, and applying quasi-elastic (QE) kinematics assuming the known incoming 
neutrino direction. Energy reconstruction in the f3B energy range is remarkably powerful, 
and the contamination of non-QE events very small, as shown in Fig. O As pointed out in 
Ref. [39], it is necessary to use a migration matrix for the neutrino energy reconstruction 
to properly handle Fermi motion smearing and the non-QE event contamination. We use 
100 MeV bins for the reconstructed energy and 40 MeV bins for the true neutrino energy. 
Four migration matrices (for z/g, z/g, t'^, z/^) are applied to signal events as well as backgrounds. 
As suggested from Fig. [2] the results using migration matrices are very similar to a Gaussian 
energy resolution. 

Backgrounds from charged pions and atmospheric neutrinos are computed with the iden- 
tical analysis chain as signal events. Charged pions generated in NC events (or in NC-like 
events where the leading electron goes undetected) are the main source of background for 
the experiment. To compute this background inclusive NC and CC events have been gen- 
erated with the f3B spectrum. Events have been selected where the only visible track is 
a charged pion above the Cerenkov threshold. Particle identification efficiencies have been 
applied to those particles. The probability for a pion to survive in water until its decay 
has been computed with Geant 3.21 and cross-checked with a Fluka 2003 simulation. This 
probability is different for positive and negative pions, the latter having a higher probability 
to be absorbed before decaying. The surviving events are background, and the reconstructed 
neutrino energy is computed misidentifying these pions as muons. Event rates are reported 
in Tab. [3l From these numbers it becomes evident that requiring the detection of the Michel 
electron provides an efficient cut to eliminate the pion background. These background rates 
are significantly smaller than quoted in Ref. [37], where pion decays were computed with 
the same probabilities as for muons and they are slightly different from those quoted in 
Ref. [48], where an older version of Nuance had been used. The numbers of Tab. [3] have 
been cross-checked by comparing the Nuance and Neugen [49] event generators, finding a 
fair agreement in background rates and energy shape. 

Also atmospheric neutrinos can constitute an important source of background [22,39,41, 
48]. This background can be suppressed only by keeping a very short duty cycle (2.2 • 10^^ 
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^He 



Ua CC 7r+ TT- 



z7„CC 7r+ TT- 



Generated ev. 
Particle ID 

Decay 



115367 557 341 
95717 204 100 
61347 107 8 



101899 674 400 
85285 240 118 
69242 120 8 



Table 3: Events for the /3B in a 4400 kt yr exposure. Vp, iy^) CC events are computed assuming 
full oscillations (P,y^_^^^ = 1), and pion backgrounds are computed from Ve{ye) CC+NC events. In the 
rows we give the number events generated within the fiducial volume ("Generated ev."), after muon particle 
identification ( "Particle ID" ) , and after applying a further identification requiring the detection of the Michel 
electron ("Decay"). 

is the target value for the /3B design study), and this in turn is one of the most challeng- 
ing bounds on the design of the Beta Beam complex. Following Ref. [48] we include the 
atmospheric neutrino background based on a Monte Carlo simulation using Nuance. Events 
are reconstructed as if they were signal neutrino events. We estimate that 5 events/year 
would survive the analysis chain in a full solar year (the /3B should run for about 1/3 of this 
period) and include these events as backgrounds in the analysis. Under these circumstances, 
the present value of the /?B duty cycle seems to be an overkill, it could be reduced by a 
factor 5 at least, see also Ref. [48] for a discussion of the effect of a higher duty cycle. 

3.3 The 3.5-GeV SPL Super Beam 

In the recent Conceptual Design Report 2 (CDR2) the foreseen Super Proton Linac (SPL) [20] 
will provide the protons for the muon production in the context of a Neutrino Factory, and 
at a first stage will feed protons to a fixed target experiment to produce an intense conven- 
tional neutrino beam ("Super Beam"). The parameters of the beam line take into account 
the optimization [25] of the beam energy as well as the secondary particle focusing and decay 
to search for z/^ — >• z/g and z/^ — >• Ve appearance as well as v^, disappearance in a mega 
ton scale water Cerenkov detector. In particular, a full simulation of the beam line from the 
proton on target interaction up to the secondary particle decay tunnel has been performed. 
The proton on a liquid mercury target (30 cm long, 7.5 mm radius, 13.546 density) has been 
simulated with FLUKA 2002.4 [50] while the horn focusing system and the decay tunnel 
simulation has been preformed with GEANT 3.21 [51]|l| 




Since the optimization requirements for a Neutrino Factory are rather different than for a 
Super Beam the new SPL configuration has a significant impact on the physics performance 
(see Ref. [25] for a detailed discussion). The SPL fiuxes of the four neutrino species (z/^, z^e, 
z/^, Ue) for the positive (z/^ beam) and the negative focusing (z/^ beam) are show in Fig. [3l The 
total number of z/^ {uy) in positive (negative) focusing is about 1.18 (0.97) x 10^^m~^y~^ with 
an average energy of 300 MeV. The Uf. (z/g) contamination in the z/^ (z/^) beam is around 
0.7% (6.0%). Following Ref. [53], the 7r° background is reduced using a tighter PID cut 

^Although there are differences between the predicted pion and kaon productions as a function of pro- 
ton kinetic energy with FLUKA 2002.4 and 2005.6, the results are consistent for the relevant energy of 
3.5 GeV. We emphasize that the pion and the kaon production cross-sections are waiting for experimental 
confirmation [52] and a new optimization would be required if their is a disagreement with the present 
knowledge. 
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Figure 3: Neutrino fluxes, at 130 km from the target with the horns focusing the positive particles (top 
panel) or the negative particles (bottom panel). The fluxes are computed for a SPL proton beam of 3.5 GeV 
(4 MW), a decay tunnel with a length of 40 m and a radius of 2 m. 



compared to the standard Super-Kamiokande analysis used in K2K, but the cuts are looser 
than for T2K. Indeed, at SPL energies the 7i° background is less severe than for T2HK. This 
is because the resonant cross section is suppressed, and the produced pions have an energy 
where the angle between the two gammas is very wide, leading to a small probability that 
the two gamma rings overlap. This results in a higher signal efficiency of SPL compared 
to T2HK (70% against 40%) and a smaller rate of ir" background. The Michel electron is 
required for the /i identification. For the Vf. channel the background consists roughly 

of 90% z/e Vf. CC interactions, 6% 7r° from NC interactions, 3% miss identified muons from 
^ii ~^ CC, and 1% z/g v^. CC interactions. For the Vf. channel the contributions 

to the background are 45% v^. Vf, CC interactions, 35% CC interactions, 18% 

7r° from NC interactions and 2% miss identified muons from z/^ — > CC. In addition we 
include the events from the contamination of "wrong sign" muon-neutrinos due to —>■ z/g 
(z/^ — >■ z/g) oscillations in the neutrino (antineutrino) mode. We have checked that with the 
envisaged duty cycle of 2.4 x 10~^ the background from atmospheric neutrinos is negligible 
for the SPL. 

Considering the signal over square-root of background ratio, the 3.5 GeV beam energy is 
more favorable than the original 2.2 GeV option. Compared to the fiuxes used in Refs. [38,53] 
the gain is at least a factor 2.5 and this justifies to reconsider in detail the physics potential 
of the SPL Super Beam. Both the appearance and the disappearance channels are used. For 
the spectral analysis we use 10 bins of 100 MeV in the interval < < 1 GeV, applying 
the same migration matrices as for the /3B to take into account properly the neutrino energy 
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Figure 4: Comparison of the fluxes from SPL and /3B. 

reconstruction. As ultimate goal suggested in Ref. [16] a 2% systematical error is used 
as default both for signal and background, this would be achieved by a special care of the 
design of the close position. However, we discuss also how a 5% systematical error affects the 
sensitivities. Using neutrino cross- sect ions on water from Ref. [33], the number of expected 
charged current is about 98 per kt yr. In Fig. H] we compare the fluxes from the SPL to 
the one from the /3B. 

3.4 The atmospheric neutrino analysis 

The simulation of atmospheric neutrino data in MEMPHYS is based on the analysis pre- 
sented in Ref. [30], with the following differences: 

• We replace the neutrino fluxes at Kamioka with those at Gran Sasso. We use the 
Honda calculations [54], which unfortunately are not yet available for the Frejus site. 
However, since the fluxes increase with the geomagnetic latitude and Frejus is northern 
than Gran Sasso, our choice is conservative. 

• We take into account the specific geometry of the MEMPHYS detector. This is partic- 
ularly important to properly separate fully contained from partially contained events, 
as well as stopping muon from through-going muon events. 

• We divide our total data sample into 420 different bins: fully contained single-ring 
events, further subdivided according to flavor (e-like or /x-like), lepton momentum (8 
bins: 0.1-0.3, 0.3-0.5, 0.5-1, 1-2, 2-3, 3-5, 5-8, 8-oo GeV) and lepton direction (20 
bins in zenith angle); fully contained multi-ring events, further subdivided according 
to flavor (e-like or /i-like), reconstructed neutrino energy (3 bins: 0-1.33, 1.33-5, 5- 
oo GeV) and lepton direction (10 bins in zenith angle); partially contained /z-like events, 
divided into 20 zenith bins; up-going muons, divided into stopping and through-going 
events, and in 10 zenith bins each. 
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• We include in our calculations the neutral-current contamination of each bin. To 
this extent we assume that the ratio between neural-current and unoscillated charged- 
current events in MEMPHYS is the same as in Super-Kamiokande, and we take this 
ratio from Ref. [3]. 

• We consider also multi-ring events, which we define as fully contained charged-current 
events which are not tagged as single-ring. Again, we assume that the survival efficiency 
and the NC contamination are the same as for Super-Kamiokande [3]. 



The expected number of contained events is given by: 



I'oo I'+l POO pEi, r~^^ /'^''^ 

7V6(cJ) = ATNC + ntstT dh dc, dE, / dEi / dc„ / d^a 

dfdcdh ^^"' ^ I '^^ ''"^ '^"^^ ' 

where P^^jj {P^^jj) is the Va — ^ ^0 {^a conversion probability for given values of 

the neutrino energy Ey^ the cosine of the angle between the incoming neutrino and the 
vertical direction, the production altitude /i, and the neutrino oscillation parameters uj. 
We calculate the conversion probability numerically in the general three-flavor framework 
taking into account matter effects from a realistic Earth density proflle. Further, Nf^ is the 
neutral-current background for the bin 6, ntgt is the number of targets, T is the experiment 
running time, {^a) ^^e flux of atmospheric neutrinos (antincutrinos) of type a, and 
a~p (ap) is the charged-current neutrino- (antineutrino-) nucleon interaction cross section. 
The variable Ei is the energy of the final lepton of type P, while Ca and if a parametrize the 
opening angle between the incoming neutrino and the final lepton directions as determined 
by the kinematics of the neutrino interaction. Finally, gives the probability that a charged 
lepton of type /3, energy Ei and direction ci contributes to the bin b. 

Up-going muon events are calculated as follows: 



r(X) p+l POO rEv pE'^ p+1 p2'K 

Nb{(^) = ProckT J] dh dc, / dE, / dE"^ / rfEj" / dCa / d<fa 

a,± J —1 J Ejnin J Ejniji J Ejnin J —1 JO 

X El") AUEl:\ c,(c,, c„, ,;.)) , (3) 

where Prock is the density of targets in standard rock, i?rock is the effective muon range [55] 
for a muon which is produced with energy and reaches the detector with energy E^^, and 
is the effective area for the bin b. The other variables and physical quantities are the 
same as for contained events. 

The statistical analysis is based on the pull method, as described in Ref. [56]. In our 
analysis we include three different kind of experimental uncertainties: Flux uncertainties: 
total normahzation (20%), tilt factor (5%), zenith angle (5%), v/V ratio (5%), and ji/e ratio 
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(5%); cross-section uncertainties: total normalization (15%) and /i/e ratio (1%) for each 
type of charged-current interaction (quasi-elastic, one-pion production, and deep-inelastic 
scattering), and total normalization (15%) for the neutral-current contributions; systematic 
uncertainties: same as in previous analyses, details are given in the Appendix of Ref. [56]. 
In addition, we assume independent normalization uncertainties (20%) for e-like and /i-like 
multi-ring events. Since we are dividing our data sample into a large number of bins, it is 
important to use Poisson statistics as some of the bins contain only a few number of events. 
We therefore write our as: 



X^ioj) = min 



Tyex 



(4) 



where the number of events for a given value of the pulls ^ is given by: 

Nl^uj, = Nl\u;) exp ^ <(^) . (5) 

The use of an exponential dependence on the pulls in Eq. (JSj), rather than the usual linear 
dependence, ensures that the theoretical predictions remain positive for any value of the 
pulls, thus avoiding numerical inconsistencies during the pull minimization procedure. 



4 Degeneracies 

A characteristic feature in the analysis of future LBL experiments is the presence of parameter 
degeneracies. Due to the inherent three-flavor structure of the oscillation probabilities, for a 
given experiment in general several disconnected regions in the multi-dimensional space of 
oscillation parameters will be present. Traditionally these degeneracies are referred to in the 
following way: 

• The intrinsic or (5cp, ^^is) -degeneracy [57]: For a measurement based on the — »• 
z/e oscillation probability for neutrinos and antineutrinos two disconnected solutions 
appear in the (^cp, ^is) plane. 

• The hierarchy or sign(Am3]^)-degeneracy [58]: The two solutions corresponding to the 
two signs of Amg]^ appear in general at different values of 6cp and ^13. 

• The octant or 6'23-degeneracy [59]: Since LBL experiments are sensitive mainly to 
sin^ 26^23 it is difficult to distinguish the two octants 6*23 < 7r/4 and 6'23 > vr/4. Again, 
the solutions corresponding to 6'23 and 7r/2 — 6*23 appear in general at different values 
of Sep and ^13. 

This leads to an eight-fold ambiguity in 6*13 and 6cp [60], and hence degeneracies provide a 
serious limitation for the determination of di^, 6cp, and the sign of Arrig^. Recent discussions 
of degeneracies can be found for example in Refs. [30,32,34,61]; degeneracies in the context 
of CERN-Frejus jSB and SPL have been considered previously in Ref. [38]. In Fig. Owe 
illustrate the effect of degeneracies for the /3B, SPL, and T2HK experiments. Assuming the 
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Figure 5: Allowed regions in sin^ 2^13 and ^cp for LBL data alone (contour lines) and LBL+ATM data 
combined (colored regions). H*''/^''(0*''/'^'') refers to solutions with the true/wrong mass hierarchy (octant 
of ^23)- The true parameter values are Sep — — O.SSvr, sin^ 2(?i3 = 0.03, sin^ 623 = 0.6, and the values from 
Eq. IT} for the other parameters. The running time is {bv + 5P) yrs for /3B and {2^ + 89) yrs for the Super 
Beams. 
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Figure 6: Resolving degeneracies in SPL by successively using the appearance rate measurement, dis- 
appearance channel rate and spectrum, spectral information in the appearance channel, and atmospheric 
neutrinos. Allowed regions in sin^ 26*13 and 6cp are shown at 95% CL, and H*''/™(0*'/™'') refers to solu- 
tions with the true/wrong mass hierarchy (octant of 023)- The true parameter values are ^cp = — 0-857r, 
sin^ 2^13 — 0.03, sin^ ^23 = 0.6, and the values from Eq. ^ for the other parameters. The running time is 
(2i^ + 8D) yrs. 



true parameter values 5cp = — O.SSvr, sin^26'i3 = 0.03, sin^ 6*23 = 0.6 we show the allowed 
regions in the plane of sin^ 29i3 and 6cp taking into account the solutions with the wrong 
hierarchy and the wrong octant of 6*23. 

As visible in Fig. [5] for the Super Beam experiments SPL and T2HK there is only a four- 
fold degeneracy related to sign(Am3]^) and the octant of 623, whereas the intrinsic degeneracy 
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can be resolved. Several pieces of information contribute to this effect, as we illustrate at 
the example of SPL in Fig. [HI The dashed curves in the left panel of this figure show 
the allowed regions for only the appearance measurement (for neutrinos and antineutrinos) 
without spectral information, i.e., just a counting experiment. In this case the eight-fold 
degeneracy is present in its full beauty, and one finds two solutions (corresponding to the 
intrinsic degeneracy) for each choice of sign(Am|^) and the octant of ^23- Moreover, the 
allowed regions are relatively large. For the thin solid curves the information from the 
disappearance rate is added. The main effect is to decrease the size of the allowed regions 
in sin^2^i3. This is especially pronounced for the solutions involving the wrong octant of 
6*23, since these solutions are strongly affected by an uncertainty in 6^23 which gets reduced 
by the disappearance information. Using in addition to the disappearance rate also the 
spectrum again decreases the size of the allowed regions, however, still all eight solutions are 
present (compare dashed curves in the right panel). The most relevant effect comes from the 
inclusion of spectral information in the appearance channel, as visible from the comparison 
of the dashed and thick-solid curves in the right panel of Fig. [6l The intrinsic degeneracy gets 
resolved and only four solutions corresponding to the sign and octant degeneracies are left 
see, e.g., Refs. [30,32,34])! Note that the thick curves in the right panel of Fig. [6] correspond 
to the regions show in Fig. for the SPL. Finally, by the inclusion of information from 
atmospheric neutrinos all degeneracies can be resolved in this example, and the true solution 
is identified at 95% CL (see Sec. 16.21 and Ref. [30] for further discussions of atmospheric 
neutrinos) . 

Concerning the jSB one observes from Fig. [5] that in this case the (5cp; 6'i3)-degeneracy 
cannot be resolved and one has to deal with eight distinct solutions. One reason for this 
is the absence of precise information on |Am|^| and sin^ 26'23 which is provided by the z/^ 
disappearance in Super Beam experiments but is not available from the /5B. If external 
information on these parameters at the level of 3% is included the allowed regions in Fig. [5] 
are significantly reduced. However, still all eight solutions are present, which indicates that 
for the /?B spectral information is not efficient enough to resolve the (5cp; ^i3)-degeneracy, 
and in this case only the inclusion of atmospheric neutrino data allows a nearly complete 
resolution of the degeneracies. 

An important observation from Fig. is that degeneracies have only a very small impact 
on the CP violation discovery, in the sense that if the true solution is CP violating also 
the fake solutions are located at CP violating values of 6cp- Indeed, since for the relatively 
short baselines in the experiments under consideration matter effects are very small, the 
sign(Am|;^)-degenerate solution is located within good approximation at (5^p ^ tt — 6cp [58]. 
Therefore, although degeneracies strongly affect the determination of 613 and 6cp they have 
only a small impact on the CP violation discovery potential. Furthermore, as clear from 
Fig. [5] the sign(Am3^) degeneracy has practically no effect on the ^13 measurement, whereas 
the octant degeneracy has very little impact on the determination of 5cp- 

Fig. [5] shows also that the fake solutions occur at similar locations in the (sin^ 26'i3, <^cp) 
plane for /3B and SPL. Therefore, as noted in Ref. [38], in this sense the two experiments are 
not complementary, and the combination of 10 years of (3B and SPL data is not very effective 
in resolving degeneracies. This is obvious since the baseline is the same and the neutrino 

^The inclusion of spectral information might be the source of possible differences to previous studies, see 
e.g. Ref. [38]. 
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95% CL regions for the (H O ), (H O ), (H O ), iH O ) solutions 
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Figure 7: Allowed regions in sin^ 26*13 and <^cp for 5 years data (neutrinos only) from /?B, SPL, and the 
combination. h*''/'"''(0*''/^'') refers to solutions with the true/wrong mass hierarchy (octant of 6*23). For the 
colored regions in the left panel also 5 years of atmospheric data are included; the solution with the wrong 
hierarchy has Ax^ = 4. The true parameter values are (5cp — — O.SStt, sin^ 26*13 — 0.03, sin^ 6*23 = 0.6, 
and the values from Eq. ^ for the other parameters. For the (JQ only analysis (middle panel) an external 
accuracy of 2% (3%) for |Ar7i|i| (^^23) has been assumed, whereas for the left and right panel the default 
value of 10% has been used. 

energies are similar. Note however, that the /3B looks for Vf. —>■ appearance, whereas in 
SPL the T-conjugate channel z/^ — > is observed. Assuming CPT invariance the relation 
^uij — Pp^^p^ holds, which implies that the antineutrino measurement can be replaced by 
a measurement in the T-conjugate channel. Hence, if /3B and SPL experiments are available 
simultaneously the full information can be obtained just from neutrino data, and in principle 
the (time consuming) antineutrino measurement is not necessary. As shown in Fig. [7] the 
combination of 5 yrs neutrino data from the f3B with 5 yrs of neutrino data from SPL leads 
to a result very close to the 10 yrs neutrino+antineutrino data from one experiment alone. 
Hence, if /3B and SPL experiments are available simultaneously the data taking period is 
reduced approximately by a factor of 2 with respect to a single experiment. This synergy is 
discussed later in Sec. 16.11 in the context of the ^13 and CP violation discovery potentials. 

5 Physics potential 

5.1 Sensitivity to the atmospheric parameters 

The z/^ disappearance channel available in the Super Beam experiments SPL and T2HK 
allows a precise determination of the atmospheric parameters lAmg^^l and sin^2^23; see, e.g., 
Refs. [62-64] for recent analyses). Fig. [S] illustrates the improvement on these parameters by 
Super Beam experiments with respect to the present knowledge from SK atmospheric and 
K2K data. We show the allowed regions at 99% CL for T2K-I, SPL, and T2HK, where in 
all three cases 5 years of neutrino data are assumed. T2K-I corresponds to the phase I of 
the T2K experiment with a beam power of 0.77 MW and the Super-Kamiokande detector 
as target [16]. In Tab. IHwe give the corresponding relative accuracies at 3a for lAmg^l and 
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Figure 8: Allowed regions of Am§i and sin^ 023 at 99% CL (2 d.o.f.) after 5 yrs of neutrino data 
taking for SPL, T2K phase I, T2HK, and the combination of SPL with 5 yrs of atmospheric neutrino 
data in the MEMPHYS detector. For the true parameter values we use Am|i = 2.2 (2.6) x lO^^ eV^ and 
sin^ 6*23 = 0.5 (0.37) for the test point 1 (2), and 6*13 = and the solar parameters as given in Eq. ((T]). The 
shaded region corresponds to the 99% CL region from present SK and K2K data [12]. 





True values 


T2K-I 


SPL 


T2HK 




2.2 ■ 10-3 eV^ 


4.7% 


3.2% 


1.1% 


sin^ ^23 


0.5 


20% 


20% 


6% 


Am2, 


2.6 ■ 10"3 eV^ 


4.4% 


2.5% 


0.7% 


sin^ 6^23 


0.37 


8.9% 


3.1% 


0.8% 



Table 4: Accuracies at Scr on the atmospheric parameters |Am|]^| and sin^ 023 for 5 years of neutrino 
data from T2K-I, SPL, and T2HK for the two test points shown in Fig. [5] {O^i^'' = 0). The accuracy for a 
parameter x is defined as (a:"PP'='' - a;'°^°'')/(2a;'™'=), where x"pp" (x'°"") is the upper (lower) bound at 3cr 
for 1 d.o.f. obtained by projecting the contour Ax^ = 9 onto the x-axis. For the accuracies for test point 2 
the octant degenerate solution is neglected. 

sin^ 6'23- 

From the figure and the table it becomes evident that the T2K setups are very good in 
measuring the atmospheric parameters, and only a modest improvement is possible with SPL 
with respect to T2K phase I. T2HK provides an excellent sensitivity for these parameters, 
and for the example of the test point 2 sub-percent accuracies are obtained at 3ct. The dis- 
advantage of SPL with respect to T2HK is the limited spectral information. Because of the 
lower beam energy nuclear Fermi motion is a severe limitation for energy reconstruction in 
SPL, whereas in T2K the somewhat higher energy allows an efficient use of spectral informa- 
tion of quasi-elastic events. Indeed, due to the large number of events in the disappearance 



16 



channel (cf. Tab. [2]) the measurement is completely dominated by the spectrum, and even 
increasing the normalization uncertainty up to 100% has very little impact on the allowed 
regions. The effect of spectral information on the disappearance measurement is discussed 
in some detail in Ref. [64]. 

In the interpretation of the numbers given in Tab.Hlone should consider that at accuracies 
below 1% systematics might become important, which are not accounted for here. We do 
include the most relevant systematics (see Sees. [2] and [3]), however, at that level additional 
uncertainties related to, for example, the spectral shapes of signal and/or background, or 
the energy calibration might eventually limit the accuracy. 

For the test point 1, with maximal mixing for 62^1 rather poor accuracies of ~ 20% for 
T2K-I and SPL, and 6% for T2HK are obtained for sin^6'23. The reason is that in the 
disappearance channel sin^ 2^23 is measured with high precision, which translates to rather 
large errors for sin^ ^23 if ^23 = 7r/4 [63]. For the same reason it is difficult to resolve the 
octant degeneracy, and for the test point 2, with a non-maximal value of sin^ ^23 = 0.37, 
for all three LBL experiments the degenerate solution is present around sin^ 62Z = 0.63. As 
pointed out in Refs. [65, 66] atmospheric neutrino data may allow to distinguish between 
the two octants of ^23- If 5 years of atmospheric neutrino data in MEMPHYS are added to 
the SPL data, the degenerate solution for the test point 2 can be excluded at more than 5a 
and hence the octant degeneracy is resolved in this example, see Sec. l6.2l for a more detailed 
discussion. 

5.2 The ^^13 discovery potential 

If no finite value of 613 is discovered by the next round of experiments an important task of the 
experiments under consideration here is to push further the sensitivity to this parameter. In 
this section we address this problem, where we use to following definition of the ^13 discovery 
potential: Data are simulated for a finite true value of sin^ 2^13 and a given true value for 
(5cp. If the Ax^ of the fit to these data with ^13 = is larger than 9 the corresponding true 
value of 6'i3 "is discovered at 3a" . In other words, the 3a discovery limit as a function of the 
true 6cp is given by the true value of sin^ 26*13 for which Ax^{0i3 = 0) = 9. In the fitting 
process we minimize the Ax^ with respect to 612, 623, Amfj, and Am^i, and in general one 
has to test also for degenerate solutions in sign(Am|]^) and the octant of 6*23. 

The discovery limits are shown for /3B, SPL, and T2HK in Fig. [9l One observes that 
SPL and T2HK are rather similar in performance, whereas the /3B with our standard 
fluxes performs significantly better. For all three facilities a guaranteed discovery reach 
of sin^ 26*13 — 5 X 10~^ is obtained, irrespective of the actual value of 6cp, however, for 
certain values of 6cp the sensitivity is significantly improved. For SPL and T2IIK discovery 
limits around sin^ 26*13 — 10~^ are possible for a large fraction of all possible values of 6cp, 
whereas for our standard /3B a sensitivity below sin^ 26*13 = 4 x 10"'^ is reached for 80% of 
all possible values of 6cp- If 10 years of data from /5B and SPL are combined the discovery 
limit is dominated by the /3B. Let us stress that the /3B performance depends crucially on 
the neutrino fiux intensity, as can be seen from the dashed curves in Fig. [HI which has been 
obtained by reducing the number of ion decays/yr by a factor of two with respect to our 
standard values given in Tab. [H In this case the sensitivity decreases significantly, but still 
values slightly better than from the Super Beam experiments are reached. 
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Figure 9: 3<t discovery sensitivity to sin^ 26*13 for /3B, SPL, and T2HK as a function of the true value of 
6cp (left panel) and as a function of the fraction of all possible values of Sqp (right panel). The running 
time is (biy + 5D) yrs for /3B and (2i^ + 8D) yrs for the Super Beams. The width of the bands corresponds 
to values for the systematical errors between 2% and 5%. The black curves correspond to the combination 
of /3B and SPL with 10 yrs of total data taking each for a systematical error of 2%, and the dashed curves 
show the sensitivity of the /3B when the number of ion decays/yr are reduced by a factor of two with respect 
to the values given in Tab. [TJ 



The peak of the sensitivity curves around Sqp ~ vr appears due to the interplay of neutrino 
and antineutrino data. For the Super Beams neutrino (antineutrino) data are most sensitive 
in the region n < ^cp ^ 27r (0 < 6cp < vr), and opposite for the /3B, compare also Fig. [THin 
Sec. 16.11 The particular shape of the sensitivity curves emerges from the relative location of 
the corresponding curves for neutrino and antineutrino data, which is controlled by the L/ Ey 
value where the experiment is operated and the value of |Am3^|. The fact that the peak is 
most pronounced for the /?B follows from the somewhat smaler L/ of the /3B compared to 
the Super Beams, whereas the shapes for SPL and T2HK are similar because of the similar 
L/Ey values. 

In Fig. [H] we illustrate also the effect of systematical errors on the ^13 discovery reach. 
The lower boundary of the band for each experiment corresponds to a systematical error of 
2%, whereas the upper boundary is obtained for 5%. These errors include the (uncorrelated) 
normalization uncertainties on the signal as well as on the background, where the crucial 
uncertainty is the error on the background. We find that the /3B is basically not affected 
by these errors, since the background has a rather different spectral shape (strongly peaked 
at low energies) than the signal. The fact that T2HK is relatively strongly affected by the 
actual value of the systematics can by understood by considering the ratio of signal to the 
square- root of the background using the numbers of Tab. [2l We shall discuss this issue in 
more detail in the next section in the context of the CP violation discovery reach. 

Let us remark that the sensitivities are practically not affected by the sign(Am3^)- 
degeneracy. This is easy to understand, since the data is fitted with 6^13 = 0, and in this case 
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Figure 10: 3a discovery sensitivity to sin^ 26113 for the SPL as a function of the true value of Sep for 
sin^ ^2™° — 0-6 and true values for the other parameters as given in Eq. ((T]). The running time is (2zy+8P) yrs. 



both mass hierarchies lead to very similar event rates. If the inverted hierarchy is used as 
the true hierarchy, the peak in the discovery limit visible in the left panel of Fig. [9] around 
Sqp ~ 71 moves to ^cp ~ 0. However, the characteristic shape of the curves, and in particular, 
the sensitivity as a function of the 5cp-fraction shown in the right panel are hardly affected 
by the sign of the true Am|]^. In case of a non-maximal value of 623 the octant- degeneracy 
has a minor impact on the ^13 discovery potential, as illustrated in Fig. [TU]for the SPL. We 
show the discovery limit obtained with the true and the fake octant of 623 for a true value 
of sin^^23 = 0.6. Let us note that for true values of sin^^23 > 0.5 the octant-degenerate 
solution leads to a worse sensitivity to 6^13 (see figure), whereas for sin^6'23 < 0.5 the fake 
solution does not affect the 613 discovery, since in this case the sensitivity is increased. 

5.3 Sensitivity to CP violation 

In case a finite value of ^13 is established it is important to quantitatively assess the discovery 
potential for leptonic CP violation (CPV). The CP symmetry is violated if the complex phase 
(5cp is different from and vr. Therefore, CPV is discovered if these values for 6cp can be 
excluded. We evaluate the discovery potential for CPV in the following way: Data are 
calculated by scanning the true values of sin^ 2^13 and ^cp- Then these data are fitted with 
the CP conserving values Sep = and ^cp = ^r, where all parameters except ^cp are varied 
and the sign and octant degeneracies are taken into account. If no fit with Ax^ < 9 is 
found CP conserving values of 6cp can be excluded at 3a for the chosen values of ^^'p'' and 
sin^ 2^f3"^ 

The CPV discovery potential for /3B, SPL, and T2HK is shown in Fig. [TTl As in the case 
of the ^13 sensitivity we find that SPL and T2HK perform rather similar, whereas the /3B has 
significantly better sensitivity if our adopted numbers of ion decays per year can be achieved. 
For systematical errors of 2% maximal CPV (for ^^^p'^ = 7i/2, 37r/2) can be discovered at 
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Figure 11: CPV discovery potential for /3B, SPL, and T2HK: For parameter values inside the ellipse- 
shaped curves CP conserving values of ^cp can be excluded at 3a (Ax^ > 9). The running time is (5z^+5P) yrs 
for /3B and {2iy + Sv) yrs for the Super Beams. The width of the bands corresponds to values for the 
systematical errors from 2% to 5%. The dashed curves show the sensitivity of the /3B when the number of 
ion decays/yr are reduced by a factor of two with respect to the values given in Tab.[l]for 2% systematics. 

3(T down to sin2 20i3 ~ 8.8(6.6) x 10"^ for SPL (T2HK), and 80^2^13 ~ 2 x 10"^ for the 
(3B. This number for the (3B is increased by a factor 3 if the fluxes are reduced to half of 
our nominal values. The best sensitivity to CPV is obtained for all three facilities around 
sin2 2^13 ~ 10-2. For this value CPV can be established for 78%, 73%, 75% of all values of 
Sep for PB, SPL, T2HK, respectively (again for systematics of 2%). 

The widths of the bands in Fig. [H] corresponds to different values for systematical errors. 
The curves which give the best sensitivities are obtained for systematics of 2%, the curves 
corresponding to the worst sensitivity have been computed for systematics of 5%. We change 
the uncertainty on the signal as well as on the background, however, it turns out that the most 
relevant uncertainty is the background normalization. We find that the impact of systematics 
is very small for the jSB. The reason for this is that the spectral shape of the background 
in the f3B (from pions and atmospheric neutrinos) is very different from the signal, and 
therefore they can be disentangled by the fit of the energy spectrum. For the Super Beams 
the background spectrum is more similar to the signal, and therefore an uncertainty on the 
background normalization might have a strong impact on the sensitivity, as visible from the 
SPL and T2HK curves in Fig. [Til In particular T2HK is strongly affected, and moving from 
2% to 5% uncertainy decreases the sensitivity to maximal CPV by a factor 3. 
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Figure 12: Impact of total exposure and systematical errors on the CPV discovery potential of /3B, SPL, 
and T2HK. We show the smallest true value of sin^ 2^13 for which 6cp — 7r/2 can be distinguished from 
^cp = or 6cp = TT at 3(7 (Ax^ > 9) as a function of the exposure in kt yrs (left) and as a function of the 
systematical error on the background (Tbkgr (right). The widths of the curves in the left panel corresponds 
to values of crbkgr from 2% to 5%. The thin solid curves in the left panel corresponds to no systematical 
errors. The right plot is calculated for the standard exposure of 4400 kt yrs. No systematical error on the 
signal has been assumed. 



This interesting feature can be understood in the following way. A rough measure to 
estimate the sensitivity is given by the signal compared to the error on the background. 
The latter receives contributions from the statistical error \/B and from the systematical 
uncertainty abkgr-B, where B is the number of background events and cxbkgr is the (relative) 
systematical error. Hence the importance of the systematics can be estimated by the ratio 
of systematical and statistical errors a^kgrB /^/B = crbkgr V^- Summing the numbers for 
background events in the neutrino and antineutrino channels given in Tab. [2] one finds that 
systematical errors dominate (abkgrV^ > 1) if (Tbkgr ^ 6%, 3%, 2% for /5B, SPL, T2HK, 
respectively. In the right panel Fig. [T2] we show the sensitivity to maximal CPV (as defined 
in the figure caption) as a function of Ubkgr- Indeed, the worsening of the sensitivity due to 
systematics occurs roughly at the values of o"bkgr as estimated above. For a more quantitative 
understanding of these curves it is necessary to consider the number of signal and background 
events for neutrinos and antineutrinos separately, as well as to take into account spectral 
information. 

The left panel of Fig. [12] shows the sensitivity to maximal CPV as a function of the 
exposur^ for values of Ubkgr from 2% to 5%. One can observe clearly that for the standard 

•^Note that the CPV sensitivity for the /3B with reduced fluxes from Fig. [TT]is worse than the value which 
follows from Fig. [121 The reason is that in Fig.[T2|the total exposure is scaled (mass x time), i.e., signal and 
background are scaled in the same way, whereas for the dashed curve in Fig. [Tl] only the fluxes are reduced 
but backgrounds are kept constant. 
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Sensitivity to CPV at 3a for pB 
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Figure 13: Impact of degeneracies on the CPV discovery potential for the (3B. We show the sensitivity to 
CPV at 3(7 (Ax^ > 9) computed for 4 different combinations of the true values of the hierarchy (NH or IH) 
and 623 (sin^ 623 — 0.4 or 0.6). Dashed curves are computed neglecting degeneracies in the fit. The running 
time is {bv + 5D) yrs. 

exposure of 4400 kt yrs T2HK is dominated by systematics and changing abkgr from 2% to 
5% has a big impact on the sensitivity. In contrast the CERN-MEMPHYS experiments (es- 
pecially the PB) are rather stable with respect to systematics and for the standard exposure 
they are still statistics dominated. We conclude that in T2HK systematics have to be under 
very good control, whereas this issue is less important for /3B and SPL. We have checked 
explicitly that the systematical error on the signal has negligible impact on these results. 
Therefore, we have set this error to zero for calculating Fig.[T2]to highlight the importance of 
the background error. In all other calculations also the signal error is included, in particular 
also in Fig. [TTl 

Let us remark that for the T2KK configuration (with one half of the Hyper-K detector 
mass at Kamioka and the second half at the same off-axis angle in Korea) the problem of 
systematics might be less severe than for T2IIK, since both detectors observe the same flux 
and background. Note however, that an important issue for the CPV sensitivity is whether 
systematics between neutrino and antineutrino data are correlated or not. We have checked 
that the worse sensitivies for T2HK shown in Fig. [TT] compared to the results obtained in 
Refs. [34, 35] can be traced back to the fact that in these papers neutrino and antineutrino 
systematics are correlated, whereas we consider them to be independent. Note that our 
approach is conservative, and the assumption of uncorrelated errors has been adopted also 
for the CERN-MEMPHYS experiments. 
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Finally, in Fig. [13] we illustrate the impact of degeneracies, as well as the true hierarchy 
and ^23-octant on the CPV sensitivity. Curves of different colors correspond to the four 
different choices for sign(Am3j^) and the ^23-octant of the true parameters. For the solid 
curves the simulated data for each choice of true sign(Am3^) and 6'23-octant are fitted by 
taking into account all four degenerate solutions, i.e., also for the fit all four combinations of 
sign(Am|]^) and 6'23-octant are used. One observes from the figure that the true hierarchy and 
octant have a rather small impact on the jSB CPV sensitivity, in particular the sensitivity to 
maximal CPV is completely independent. The main effect of changing the true hierarchy is to 
exchange the behavior between < 5cp < 180° and 180° < 5cp < 360°. For sin2 2^i3 < 10"^ 
the sensitivity gets slightly worse if ^23^'^ > it/A compared to ^2™° < vr/4. 

The dashed curves in Fig. [13] are computed without taking into account the degeneracies, 
i.e., for each choice of true sign(Am|^) and ^23-octant the data are fitted only with this 
particular choice. The effect of the degeneracies becomes visible for large values of 61^ . 
Note that this is just the region where they can be reduced by a combined analysis with 
atmospheric neutrinos (see Sec. 16.21 or Ref. [30]). 

6 Synergies provided by the CERN-MEMPHYS facilities 

6.1 Combining Beta Beam and Super Beam 

In this section we discuss synergies which emerge if both jSB and SPL are available. The main 

difference between these two beams is the different initial neutrino flavor, 'z/'g for /?B and *z/'^ 
for SPL. This implies that at near detectors all relevant cross sections can be measured. In 
particular, the near detector of the /3B will measure the cross section for the SPL appearance 
search, and vice versa. If both experiments run with neutrinos and antineutrinos all possible 
transition probabilities are covered: Pi^^^^^, Pp^^p^, P^^^^,^, and Pp^^p^. Together with the 
fact that matter effects are very small because of the relatively short baseline, this means 
that in addition to CP also direct tests of the T and CPT symmetries are possible. 

However, if the CPT symmetry is assumed in principle all information can be obtained 
just from neutrino data because of the relations Pp^^p^ = Puf^^u^ and Pp^^p^ = Pu^^v^- 
As mentioned already in Sec. H] this implies that (time consuming) antineutrino running 
can be avoided. We illustrate this synergy in Figs. [H] and [13 In Fig. [H] we show the ^13 
discovery potential of 5 years of neutrino data from /3B and SPL. From the left panel the 
complementarity of the two experiments is obvious, since each of them is most sensitive in a 
different region of 5cp. (As expected from general properties of the oscillation probabilities 
the sensitivity curves of /?B and SPL are approximately related by the transformation 5cp ~^ 
27r — (5cp-) Combining these two data sets results in a sensitivity slightly better than from 10 
years (2i/+8i/) of T2IIK data. As visible in Fig. [15] also for the CPV discovery this synergy 
works and 5 years of neutrino data from /3B and SPL lead to a similar sensitivity as 10 years 
of T2HK. 

6.2 Resolving degeneracies with atmospheric neutrinos 

It was pointed out in Ref. [30] that for LBL experiments based on mega ton scale water 
Cerenkov detectors data from atmospheric neutrinos (ATM) provide an attractive method 
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3a discovery of a non-zero 9 within 5 yrs 




Figure 14: Discovery potential of a finite value of sin^ 2^13 at 3a {Ax^ > 9) for 5 yrs neutrino data from 
/3B, SPL, and the combination of f3B + SPL compared to 10 yrs data from T2HK (2 yrs neutrinos + 8 yrs 
antincutrinos). 
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3a discovery of CP violation: Ax (6„ = 0, 7t) = 9 
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Figure 15: Sensitivity to CPV at 3ct (Ax^ > 9) for combining 5 yrs neutrino data from /3B and SPL 
compared to 10 yrs data from T2HK (2 yrs neutrinos + 8 yrs antincutrinos). 

to resolve degeneracies. Atmospheric neutrinos are sensitive to the neutrino mass hierarchy if 
^13 is sufficiently large due to Earth matter effects, mainly in multi-GeV e-like events [67-69]. 



24 



2a sensitivity to normal hierarchy from LBL + ATM data 
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Figure 16: Sensitivity to the mass hierarchy at 2(t (Ax^ = 4) as a function of the true values of sin^ 2^13 
and 6cp (left), and the fraction of true values of Sep (right). The solid curves are the sensitivities from the 
combination of long-baseline and atmospheric neutrino data, the dashed curves correspond to long-baseline 
data only. The running time is (5;^ -I- bP) yrs for /3B and (2z/ -I- 8P) yrs for the Super Beams. For comparison 
we show in the right panel also the sensitivities of NOi^A and N0i^A-|-T2K extracted from Fig. 13.14 of 
Ref [17]. For the curve labeled "NOi^A (p.dr.)+T2K@4 MW" a proton driver has been assumed for NOi^A 
and the T2K beam has been up-graded to 4 MW, see Ref [17] for details. 

Moreover, sub-GeV e-like events provide sensitivity to the octant of 623 [65, 66, 70] due to 
oscillations with Am2i (see also Ref. [71] for a discussion of atmospheric neutrinos in the 
context of Hyper-K). Following Ref. [30] we investigate here the synergy from a combination 
of LBL data from pB and SPL with ATM data in the MEMPHYS detector. Technical details 
are given in Sec. 13. 4[ 

The effect of degeneracies in LBL data has been discussed in Sec. HI see Figs. [5] and [6l 
As discussed there, for given true parameter values the data can be fitted with the wrong 
hierarchy and/or with the wrong octant of 623- Hence, from LBL data alone the hierarchy 
and the octant cannot be determined and ambiguities exist in the determination of di^ and 
6cp- If the LBL data are combined with ATM data only the colored regions in Fig. [S] survive, 
i.e., in this particular example for SPL and T2HK the degeneracies are completely lifted at 
95% CL, the mass hierarchy and the octant of 623 can be identified, and the ambiguities in 
613 and ^cp are resolved. For the /3B an island corresponding to the wrong hierarchy does 
survive at the 95% CL for 2 dof. Still, the solution with the wrong sign of Amg^^ is disfavored 
with Ax^ = 5.1 with respect to the true solution, which corresponds to 2. lex for 1 dof. Let us 
note that in Fig. Elwe have chosen a favorable value of sin^ 623 = 0.6; for values sin^ 623 < 0.5 
in general the sensitivity of ATM data is weaker [30]. 

In Fig. [16] we show how the combination of ATM-I-LBL data leads to a non-trivial sensi- 
tivity to the neutrino mass hierarchy, i.e. to the sign of Amg^^. For LBL data alone (dashed 
curves) there is practically no sensitivity for the CERN-MEMPHYS experiments (because 
of the very small matter effects due to the relatively short baseline), and the sensitivity 
of T2HK depends strongly on the true value of 6cp- However, by including data from at- 
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mospheric neutrinos (solid curves) the mass hierarchy can be identified at 2cr CL provided 
sin^ 2^13 > 0.02 — 0.03. As an example we have chosen in that figure a true value of = 7r/4. 
Generically the hierarchy sensitivity increases with increasing 623, see Ref. [30] for a detailed 
discussion. 

The sensitivity to the neutrino mass hierarchy shown in Fig. [16] is significantly improved 
with respect to our previous results obtained in Ref. [30]. There are two main reasons for 
this improved performance: First, we use now much more bins in charged lepton energy 
for fully contained single-ring event^ (compare Sec. 13. 4p . and second, we implemented also 
information from multi-ring events. This latter point is important since the relative contri- 
bution of neutrinos and antineutrinos is different for single- and multi-ring events. Therefore, 
combining single- and multi-ring data allows to obtain a discrimination between neutrino 
and antineutrino events on a statistical basis. This in turn contains crucial information on 
the hierarchy, since the matter enhancement is visible either in neutrinos or antineutrinos, 
depending on the hierarchy. 

Although /?B and SPL alone have no sensitivity to the hierarchy at all, we find that the 
combination of them does provide rather good sensitivity even without atmospheric data. 
The reason for this interesting effect is the following. Because of the rather short baseline 
the matter effect is too small to distinguish between NH and IH given only neutrino and 
antineutrino information in one channel. However, the tiny matter effect suffices to move the 
hierarchy degenerate solution to slightly different locations in the (sin^2^^i3, 6cp) plane for 

the (/?B) and e (SPL) channels (compare Fig. EI). Hence, if all four CP and 

T conjugate channels are available (as it is the case for the /9B+SPL combination) already 
the small matter effect picked up along the 130 km CERN-MEMPHYS distance provides 
sensitivity to the mass hierarchy for sin^ 2^^i3 > 0.03, or sin^ 26'i3 ^ 0.015 if also atmospheric 
neutrino data is included. 

For comparison we show in the right panel of Fig. [16] also the sensitivity of the NOi/A [17] 
experiment, and of N0z/A+T2K, where in the second case a beam upgrade by a proton driver 
has been assumed for NOz/A, and for T2K the Super-Kamiokande detector has been used 
but the beam intensity has been increased by assuming 4 MW power. More details on these 
sensitivities can be found in Ref. [17]. Let us note that in general LBL experiments with two 
detectors and/or very long baselines (> 1000 km) are a competitive method to atmospheric 
neutrinos for the hierarchy determination, see, e.g., Refs. [34-36,73,74] for recent analyses. 
In particular, in case of the T2KK extension of the T2HK experiment, the very long baseline 
to Korea allows for a determination of the mass hierarchy down to sin^ 2^13 > 0.02 {2a) 
without using atmospheric neutrinos [35]. We mention also the possibility to determine the 
neutrino mass hierarchy by using neutrino events from a galactic Super Nova explosion in 
mega ton Cerenkov detectors such as MEMPHYS, see, e.g., Ref. [75]. 

Fig. [T7I shows the potential of ATM-I-LBL data to exclude the octant degenerate solution. 
Since this effect is based mainly on oscillations with ^m\^ there is very good sensitivity even 
for 6*13 = 0; a finite value of ^13 in general improves the sensitivity [30]. From the figure 
one can read off that atmospheric data alone can can resolve the correct octant at 3(T if 
I sin^ ^23 — 0.5| > 0.085. If atmospheric data is combined with the LBL data from SPL 
or T2HK there is sensitivity to the octant for | sin^ 6'23 — 0.5| > 0.05. The improvement 

^The impact of energy binning on the hierarchy determination with atmospheric neutrinos has been 
discussed recently in Ref. [72] in the context of magnetized iron detectors. 
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Figure 17: Ax^ of the solution with the wrong octant of 823 as a function of the true value of sin^ 023. 
We have assumed a true value of 613 ~ 0, and the running time is {2i/ + 81?) yrs. 

of the octant sensitivity with respect to previous analyses [30, 66] follows from changes in 
the analysis of sub-GeV atmospheric events, where now three bins in lepton momentum are 
used instead of one. Note that since in Fig. [T7] we have assumed a true value of 6^13 = 0, 
combining the (3B with ATM does not improve the sensitivity with respect to atmospheric 
data alone. We note that the T2KK configuration provides also some sensitivity to the 
octant of 623 without referring to atmospheric neutrinos. In Ref. [35] it was found that for 
I sin^^23 ~ 0.5| > 0.12 the octant can be identified in T2KK at 3a. 

7 Summary 

In this work we have studied the physics potential of the CERN-MEMPHYS neutrino os- 
cillation project. We consider a Beta Beam (/3B) with 7 = 100 for the stored ions, where 
existing facilities at CERN can be used optimally, and a Super Beam based on an optimized 
Super Proton Linac (SPL) with a beam energy of 3.5 GeV and 4 MW power. As target we 
assume the MEMPHYS detector, a 440 kt water Cerenkov detector at Frejus, at a distance of 
130 km from CERN. The main characteristics of the experiments are summarized in Tab. [H 
The adopted neutrino fluxes are based on realistic calculations of ion production and storage 
for the f3B, and a full simulation of the beam line (particle production and decay of secon- 
daries) for SPL. Special care has be given to the issue of backgrounds, which we include by 
means of detailed event simulations and applying Super-Kamiokande particle identification 
algorithms. 

The physics potential of the /3B and SPL experiments in terms of 613 discovery reach and 
sensitivity to CP violation has been addressed where parameter degeneracies are fully taken 
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into account. The main results on these performance indicators are summarized in Figs. [9] 
and[TTl We obtain a guaranteed discovery reach of sin^ 2^13 ^ 5 x 10^^ at 3a, irrespective 
of the actual value of 6cp- For certain values of 6cp the sensitivity is significantly improved, 
and for /3B (SPL) discovery limits around sin^ 26'i3 ~ 3 (10) x 10~^ are possible for a large 
fraction of all possible values of 6cp- Maximal CP violation (for S^p*^ = 7r/2, 37r/2) can be 
discovered at 3a down to sin 26*13 ^ 2 (9) X 10"^ for f3B (SPL), whereas the best sensitivity 
to CP violation is obtained for sin^ 26*13 ~ 10~^: For sin^ 26*13 = 10~^ CP violation can 
be established at 3a for 78% (73%) of all possible true values of 6cp for /3B (SPL). We 
stress that the /3B performance in general depends crucially on the number of ion decays 
per year. The impact of the value of systematical uncertainties on signal and background 
on our results is discussed. The /3B and SPL sensitivities are compared to the ones of the 
phase II of the T2K experiment in Japan (T2HK), which is a competing proposal of similar 
size and timescale. In general we obtain rather similar sensitivities for T2HK and SPL, and 
hence the CERN-MEMPHYS experiments provide a viable alternative to T2HK. We find 
that /3B and SPL are less sensitive to systematical errors, whereas the sensitivity of T2HK 
crucially depends on the systematical error on the backgroundjf] 

Assuming that both (3B and SPL experiments are available, we point out that one can 
benefit from the different oscillation channels Ue z/^ for /5B and for SPL, since 

by the combination of these channels the time intensive antineutrino measurements can be 
avoided. We show that 5 years of neutrino data from /3B and SPL lead to similar results as 2 
years of neutrino plus 8 years of antineutrino data from T2HK. Furthermore, we discuss the 
use of atmospheric neutrinos in the MEMPHYS detector to resolve parameter degeneracies 
in the long-baseline data. This effect leads to a sensitivity to the neutrino mass hierarchy 
at 2a CL for sin^2^i3 > 0.025 for /3B and SPL, although these experiments alone (without 
atmospheric data) have no sensitivity at all. The optimal hierarchy sensitivity is obtained 
from combining /?B+SPL+atmospheric data. Furthermore, the combination of atmospheric 
data with a Super Beam provides a possibility to determine the octant of ^23- 

To conclude, we have shown that the CERN-MEMPHYS neutrino oscillation project 
based on a Beta Beam and/or a Super Beam plus a mega ton scale water Cerenkov detector 
offers interesting and competitive physics possibilities and is worth to be considered as a se- 
rious option in the worldwide process of identifying future high precision neutrino oscillation 
facilities [76]. 
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